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ABSTRACT 


The  objective  oi  this  report  is  to  present  the  effect 
of  a  tuned  damper  on  a  single  degree-of-f reedom  system 
which  is  subjected  to  white  noise  excitation.  The  tuned 
damper  itself  consists  of  a  mass  connected  to  a  visco¬ 
elastic  link  which,  in  turn,  is  connected  to  the  primary 
system  under  consideration.  The  criterion  used  for  tuning 
the  damper  is  the  minimization  of  the  mean  square  response 
of  the  primary  system.  The  tuned  damper  obtained  by  use 
of  this  criterion  is  compared  to  that  obtained  from  another 
criterion  requiring  the  peaks  of  the  absolute  value  of  the 
frequency  response  function  to  be  of  equal  height. 
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1.  INTRODUCTION 


Reduction  of  the  amplitude  of  a  vibrating  system  through 
the  use  of  various  damping  devices  has  received  extensive 
study  (see  References  1,  4,  7,  8,  and  9).  A  method  of  damping 
that  is  presently  receiving  some  attention  is  that  of  using 
viscoelastic  damper  units  (see  References  3,  4,  5,  and  10). 

This  unit  consists  of  a  mass  attached  to  a  viscoelastic  link 
which  is,  in  turn,  attached  to  the  primary  system  under  con¬ 
sideration. 

Consider  a  primary  system  consisting  of  a  spring  with  a 
spring  rate  k^  which  suspends  a  mass  m^  as  shown  in  Figure  1. 

If  a  viscoelastic  damper  unit,  idealized  as  a  small  mass 
connected  to  a  spring  with  a  complex  modulus,  is  attached  to 
the  primary  system,  the  absolute  value  of  the  frequency 
response  function  exhibits  two  peaks.  The  frequency  response 
function  is  the  complex  ratio  of  the  steady-state  response 
of  the  system  to  a  sinusoidal  input  and  is  a  function  of  the 
excitation  frequency.  It  has  been  suggested  (see  References  10, 
11)  that  an  optimum  damper  would  damp  the  motion  of  the 
primary  system  such  that  the  two  peaks  of  the  absolute  value 
of  the  frequency  response  function  would  be  of  equal  height. 

This  report  outlines  an  investigation  into  another 
possible  criterion  for  optimization  of  viscoelastic  dampers. 

The  same  primary  spring-mass  system,  to  which  is  attached  a 
viscoelastic  damper  unit,  is  considered.  The  system  is 
excited  by  a  special  type  of  random  excitation,  a  white  noise. 
The  criterion  that  is  proposed  for  optimization  of  this  system 
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is  that  the  mean  square  value  of  the  weakly  stationary  random 
response  be  a  minimum.  Although  a  white  noise  is  considered, 
the  results  are  applicable  to  any  excitation  which  exhibits  a 
broad  band  flat  spectral  density  over  a  range  of  frequencies 
of  practical  interest.  See  Reference  12. 

The  method  of  residues  for  complex  variables  is  employed 
in  determining  the  mean  square  value  of  the  weakly  stationary 
response.  Numerical  results  are  presented  for  various  combi¬ 
nations  of  the  parameters  involved  and  are  compared  with  those 
corresponding  to  the  criterion  of  equal  peaks  for  the  absolute 
value  of  the  frequency  response  function.  The  numerical 
computations  were  carried  out  on  a  high  speed  computer. 
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2.  ANALYSIS 


As  is  stated  in  the  introduction,  the  criterion  that  is 
proposed  for  optimization  of  the  viscoelastic  damper  unit  is 
to  obtain  a  minimum  mean  square  value  of  the  weakly  stationary 
response  when  the  system  is  excited  by  a  white  noise.  When 
this  condition  is  satisfied,  the  damper  unit  will  be  referred 
to  as  tuned.  Since  the  input  is  random,  the  response  is  also 
random.  The  springs  and  masses  making  up  the  primary  system 
and  the  damper  unit  are  considered  to  be  deterministic. 

If  the  system  is  excited  by  a  weakly  stationary  random 
excitation,  such  as  a  white  noise,  the  response  of  the  system 
also  becomes  weakly  stationary  after  a  sufficient  amount  of 
time  has  passed  so  that  the  transient  motion  has  died  out. 

In  this  manner,  the  weakly  stationary  response  is  analogous 
to  the  steady  state  response  in  deterministic  vibration  theory. 
It  may  be  called  the  steady  state  in  the  probabilistic  sense. 
The  time  required  for  this  condition  to  occur  is  dependent 
upon  the  amount  of  damping  in  the  system..  The  greater  the 
damping,  the  sooner  the  response  becomes  weakly  stationary. 

Two  cases  involving  random  excitation  are  considered. 
Figure  1  corresponds  to  the  case  of  random  excitation  applied 
at  mass  of  the  system.  Figure  2  corresponds  to  the  case 
where  the  foundation  is  moved  in  a  random  fashion.  The  tuned 
dampers  corresponding  to  these  cases  are  compared  to  the 
tuned  dampers  corresponding  to  the  criterion  for  equal  peaks 
of  the  absolute  value  of  the  respective  frequency  response 
functions. 
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When  the  random  input  is  weakly  stationary,  the  autocor¬ 
relation  function  of  the  input  is  a  function  of  the  difference 
in  the  parametric  values 


-  ^FF ^  ^  1  -  2^ 


(1) 


The  autocorrelation  function  can  bo  expressed  as  the  founer 
transform  of  the  spectral  density;  that  is 


RFF  *  ^  1  ”  ^  2  ^  /$FF(oO)e,*'U">^l  ^  2^  du)  12 


In  order  to  compute  the  correlation  function  of  the 
response,  use  ts  made  of  the  following  relationship  between 
the  random  input  and  the  random  output 


Erx(t1)x<t2)i  = 


A 


rx2 


d  t 


1 


d  t2  EtF(t1)F(r2)l 

i  h(tt  - 1 -  r2) 


(31 


It  has  been  assumed  that  Xft)  =  X ( t )  =  0  at  t  =  0  with  prob¬ 
ability  one.  Substituting  Equations  O )  and  (2)  into  {3) 
the  correlation  function  of  the  response  can  b-^  written 

t  t  g  j~oa 

ELX<  t1)X(t2>  I  =  j  dtj  J  d  tr  2  I  du|FpfcO) 

J° 

( 4 ) 

eiuoCtj_  _t2)h(t1  -t:i>h(t2  -  t2) 


Integrating  first  on  “t  ^  and  X.  2 

ErX(t1)X(t2)J  =  (tO,tl)  X‘ltO,t2) 


e<l0fti  ”  x2*  duo 


(5) 


It 


where  the  following  notation  has  been  used: 


o 


(6) 


The  interchange  of  the  order  of  integration  in  Equation  (4) 


is  permissible  provided  that  the  function  X  (CO,t)  is 


uniformly  bounded  in  u)  .  This  condition  is  always  satisfied 
for  systems  with  positive  damping.  It  may  be  noted  that  the 
lower  limit  of  the  integral  in  Equation  (6)  can  be  extended 
to  -  oe  ,  since  h(u)  vanishes  for  negative  u.  Furthermore, 
if  the  upper  limit  t  tends  to  infinity,  then  the  right  hand 
side  of  Equation  (6)  becomes  the  frequency  response  function, 
that  is: 


£4  (60  ,  )  =  H  (U)  ) 


Since  only  the  mean  square  value  of  X(t)  in  the  weakly 
stationary  state  will  be  considered,  we  let  t^  =  tg  =  t  and 
let  t  tend  to  infinity  in  Equation  (5).  We  obtain: 


|H(ui)|  2  |  ff(60)  doO  (7) 


E[X2(t)]  = 


Equation  (7)  describes  the  relationship  between  the  weakly 
stationary  mean  square  response  of  the  system,  the  spectral 
density  of  the  random  input,  and  the  absolute  value  squared 
of  the  frequency  response  function.  The  absolute  value 
squared  of  the  frequency  response  function  prescribes  the 
fraction  of  energy  to  be  transmitted  through  the  system  at 
various  frequencies. 
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Integration  of  Equation  (7)  can  be  performed  if  the 
spectral  density  of  the  excitation  is  known.  The  random 
excitation  considered  in  this  report  is  white  noise  excitation. 
For  this  excitation,  the  spectral  density  is  a  constant.  The 
physical  interpretation  of  a  constant  spectral  density  is 
that  the  energy  content  in  the  random  forcing  function  is 
uniformly  distributed  over  the  entire  frequency  domain. 

Since  this  corresponds  to  an  infinite  mean  energy,  a  white 
noise  is  physically  impossible.  However,  if  the  absolute 
value  of  the  frequency  response  function  is  sharply  peaked 
near  the  natural  frequencies  of  the  system,  and  the  actual 
input  spectral  density  varies  slowly  in  the  neighborhood  of 
the  peaks, then  the  excitation  can  be  treated  as  white  noise 
while  computing  the  second  order  properties  of  the  response. 

The  white  noise  excitation  is  then  physically  meaningful  in 
the  sense  of  being  a  good  approximation  to  an  actual  spectral 
density  for  such  computations. 

Since  the  spectral  density  for  a  white  noise  is  a  con¬ 
stant,  K,  the  expression  for  the  mean  square  value  of  the 
weakly  stationary  response  can  be  written: 

E[X2(t)]  =  K  |H(o))  |  2  d  CO  (8) 

-  00 

The  viscoelastic  element  in  the  damper  unit  is  represented 
by  a  spring  with  a  complex  valued  stiffness, 
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kgCl  +  ),  where  is  referred  to  as  the  loss  factor.  For 

the  system  shown  in  Figure  1,  it  can  be  shown  that: 


H(CO) 


[wf  (l+i^)  -  CO4] 

1 

m, 

|w4-u^u)*(n.m)*u)TcoS  1^10^4.0*1,1 

r>il[c04,-u^tO*-(|+M)-u;^L^+  u>fc*)J  ]-  i.W|  [-u)\t01^(l+M)  +  ^U3TtOi  j 


where : 


(10) 


This  frequency  response  function  refers  to  the  displacement 
of  mass  m^.  It  is  the  response  of  this  mass  that  is  of 
primary  importance  since  this  mass  represents  the  physical 
structure . 

Equation  (8)  can  be  evaluated  by  the  method  of  residues. 
The  integrand,  which  is  a  function  of  the  real  variable  CO  , 
can  be  treated  as  a  function  of  a  complex  variable  A 
Since  the  method  of  residues  is  employed,  the  poles  of  the 
right  hand  side  of  Equation  (9)  must  be  located.  Consider 
the  following  portion  of  the  integrand  (where  CO  has  been 
replaced  by  A  ): 

tOa,  ( I  +  ^  

To  locate  the  poles,  the  denominator  is  set  equal  to  zero: 


A4-  A2[W2(1+m)+  Wi  +  t\U>2<1+M)  1  +0Ji  a^(1+  >  *  0  (12) 
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V  2 

Solving  for  A  : 

A2  =  i[U>2(l+M)  +0)2  +  l*|^2(l+M)] 

(13) 
1/2 

-  5  |  [C02(1+M)  +  U>2  +  C^C02(1+M)]2  -  4t02  602(l+i^  )  j 

This  expression  can  be  rewritten  as  follows: 

^2  ±  ^  cosl^ 

2  2  V  2  / 

(14) 

+  |[^2A  1  ^  sin(HI 


where:  D  =  \]  R2  +  I2 '  (15) 

9  =  tan  ^  ^  (16) 

R  =  A2U0^  +  Ol2602  (  2  A-4)  +  -  *^2A2  U)g  (17) 

I  =  2>^A2cO^  +^ui2t02  (^2 A-4)  (18) 

A  =  1  +  M  (19) 


Note  that  the  principal  value  of  9  is  to  be  used  in  Equation 
(14).  If  the  following  substitutions  are  made, 


E1  = 


E2  = 


F  = 
1 


F2  = 


ioj  A  +  lof 


^Tcos 

z 


(!) 

A  -  ..  -  (41 

Z  Z  \Z) 

A  +  NflT1  Sin  /  £  ^  ] 

^  [^lA  -  \TFsin  (|)] 


then  the  expression  for  A  can  be  further  simplified  to: 

A2  =  Ex  +  iT1 

A2  =  E„  +  i  F0 


(20) 

(21) 

(22) 

(23) 

(24) 

(25) 
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The  poles  of  expression  (11)  are  now  given  as  follows: 


#  * 

A  i  =  \[g^[cos  (^j-)  +  i  sin  (^2>]  (26) 

A  2  =  -fG^[cos  (^-)  +  i  sin  (\)  ]  (27) 

A  3  -  ^  Gg[cos  (^5-)  +  i  sin  (^— )  ]  (28) 

A  4  =  — Gg’lcos  (^2")  +  *■  s*n  (— 2~ ^  (29) 

where:  Gj^  =  'J  Ex  + 

0*  =  tan"1  *1/El  <31> 

G2  =  02) 

©**  =  tan’1  F2/E2  (33) 


*  ** 

Again,  only  the  principal  values  of  ©  and  ©  are  used  in 
Equations  (26)  thru  (29).  Or,  more  simply: 


Ax  =  a  +  ib 

(34) 

A2  =  -a  -  ib 

(35) 

A 3  =  u  +  iv 

(36) 

A  4  =  -u  -  iv 

(37) 

Returning  to  Equation  (9)  it  is  seen  that  the  portion 
of  the  absolute  value  squared  of  the  frequency  response 
function  that  has  not  been  considered  is  the  complex  conjugate 
of  the  portion  that  has  just  been  considered.  It  can  be  shown 
that  the  poles  of  the  portion  yet  to  be  considered  are  of  the 


form: 

*5 

=  a  -  ib 

(38) 

=  -a  +  ib 

(39) 

X7 

=  u  -  iv 

(40) 

=  -u  +  iv 

(41) 
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If,  for  example,  a,  b,  u,  and  v  are  all  positive,  one  possible 
location  of  the  eight  values  of  A  is  shown  in  Figure  3. 

We  are  now  ready  to  perforin  the  integration  shown  in 
Equation  (8).  The  limits  on  this  integral  are  minus  infinity 
and  plus  infinity.  Replace  the  right  hand  side  of  Equation 
(8)  by  a  contour  integration  and  let  the  contour  be  as  shown 
in  Figure  3.  Application  of  the  residue  theorem  yields: 

E[X2  (t)  ]  =  2  TT  Ki(R1  +  R3  +  Rg  +  Rg)  (42) 


where  Rj is  the  residue  at  A  =  A  j .  These  residues  are 


readily  evaluated  as  follows: 


R,  = 


R„  = 


co 2  +  60 2  n  -Zu)2  A2  +  A* 


m 


T<V  V<  V  As)(  Va4>---<  ArA8) 

CO  g  +  0^2  \  2  -2^2  ^3  +  ^3 

^(A3- A^f  A3-  A2)(  A3-  A4)...(  A3-  Ag> 


CO  2  +  60^  “  ^602  ^6  +  ^ 


R6  2 


m 


\(  Ag-  Ax>  (  Ag-  A2) . . .  ( A 6-  A5)  (  Ag-  A  7)  (  Ag-  Ag) 

4  .  ,  2  o  .  .  2  i  2  i  4 


R8  “  2 


^2  +  ^2  \  ~  ZtO 2^8  +  ^8 


Ag-  (  Ag-  A2)  (  Ag-  A3) . . .  (  A8-  A?) 

To  further  facilitate  later  computations,  the  following 
substitutions  are  made: 


(43) 


(44) 


(45) 


(46) 


a  +  u  =  c 
a  -  u  =  d 
b  +  v  =  g 
b  -  v  =  h 
dh  +  dg  =  p 
eg  +  ch  =  q 
d  -  hg  =  r 


rs 

— 

pq 

t 

ps 

+ 

qr 

= 

w 

dh 

- 

dg 

- 

X 

d2 

+ 

gh 

y 

gc 

- 

ch 

z 

2 

c 

+ 

gh 

= 

n 

yn 

- 

xz 

m 

(47) 


C  -  hg  =  s 

9  2 

8ab  t  +  8a  bw  = 

2  2 
8u  vk  +  8uv  m  — 


xn  +  yz  =  k 
9  9 

8a  bt  -  Sab  w  =  0“ 
8u2vm  -  8uv2k  = 


Then  RY  +  R3  +  Rg  +  Rg  can  be  written  as  follows: 

^2+ ^2  ^^2  ^1+ ^1  ^  2+  ^2  ^  ~  ^  ^  2  ^  3+  ^  3 

-£>  +icr 

+  L 04+  ’l2"  2.  6O2  Ag+  Ag  ^2+  ^  \  ~  ^  ^2  ^8+  ^ 


(48) 


8 


p  +  io~  ^  + 1  y 

The  sum  of  these  four  terms  is  a  fraction  whose 
denominator  is: 


S'2)  (49) 


(-^>+i(T)  (  £>+i0-)(-^+itf  )(^  +  i^  )  «  (p2+C2)(^‘ 

The  numerator  is  very  lengthy.  Since  the  denominator  is  a 
real  quantity,  the  numerator  must  be  a  purely  imaginary 
quantity,  otherwise  E[X2(t)]  computed  from  Equation  (8)  would 
be  complex.  If  multiplication  of  all  the  terms  in  the  numer¬ 
ator  is  carried  out,  it  is  seen  that  the  real  terms  do  cancel. 
The  numerator  so  obtained  is  given  as  follows: 

NUMERATOR  =”^r^  \  -4abC02  2  +  4a3b  ^ 


ml  v 

3  _  /o  2  _ x2  ^2  -  3.  v*  2 


-  4ab  “  4ab^02  + 


4a  b 


-  4ab 


3  _  v  2 


C  7 

do^cr^2  +  602  ^2o-^2  -  Zu)22  acr^2  +  Zco22b2cr^2 


(50) 


2— 2<rlT 2 


+  aV^2  -  6a2bV^2  +  b4cr^  2  +  ^2  CT  ^  2  + 

-  Z  Vtf  2  +  2o&  b2c r  tf2  +  a4crtf 2  -  eaVcr  X2 
2  £ 
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+  b4  O'  2  -  +  4uJv^^|Z  -  4uvJ<^  0 

-  4uv£*>2  <^(T2  +  4u3v  T2  -  4uv3^CT2  +  CO4  p  2 

+  CO2  >|^2  ?f  2  -  Z  6l>2U^  +  £^2  v2^p2  +  u^  3"^ 

~  6u2v2ffp  2  +  v4  tfp2  +  £04  tfcr2  +  cO2^2  ffcr  2 

-  £co2  u2  cr 2  +  ^602v2 so- 2  +  u4 tfcr 2  -  su2v2  tfcr  2 


.2^-2 


3»  J 


+  v4  if  cr  2  i 

Therefore :  E[X2(t)]  «  2  ”i 

(^+(T2)(^2+  r2) 


(51) 


Computations  involving  Equation  (51)  were  carried  out  on  a 
high  speed  computer. 

The  other  case  that  is  considered  is  shown  in  Figure  2. 
It  is  similar  to  the  system  shown  in  Figure  1,  but  the 
excitation  is  the  random  motion  of  the  foundation.  Since 
the  excitation  is  again  a  white  noise,  Equation  (8)  applies 
to  this  case,  except  that  the  square  of  the  absolute  value 
of  the  frequency  response  function  is  now  given  by: 


H( A  )|2  = 


[u>l  (I  +  C 

i)  ~ 

*■  m)  -udfAz+  u)f  tO§| 

1+i.w.j 

(52) 


^ .  [A4-  di  A4  (l+M  )-Uf. A*  +  dtu>l  ]  -  LM 1  [-  at  0+M>  ^u>l] 


This  frequency  response  function  refers  to  the  difference 

between  the  displacement  of  mass  m^  and  the  displacement  of 

the  foundation,  that  is,  -  Y  =  Z. 

v  4 

The  presence  of  A  in  the  numerator  of  the  integrand 
does  not  affect  the  location  of  the  poles  of  the  integrand. 
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but  it  does  affect  the  values  of  the  residues.  It  can  be 

2 

shown  that  the  expression  for  E[Z  (t)]  can  be  written  as 
follows:  o 


E [Z  (t )  ]  = 


STT  K  r 

2<e2+a-2)(^2+lf2)  '-4abw2^ 


2  -  2 


m 


4a 


t ,  ^  '■ 

be1 


-  4ab 


-  4abt0j  +  4a3b  2  -  4ab3p'J’ 


+  a 


4  _  <a  2 


4ab  tOf; 

2“  w  ^  - - 2’ 

b*^*  +  2  +  u>?-2~^2 


4^&2  .  ,.N2^2a-^2  _  <^£<o2a20.^2  +  £^^2^  2 


lc ry  -  6aVor^z  +  b^cry  +  LO*crS  -  + 


-  242^2  +  -  +  a 


2-2-  *  2  -4Q-tf  2  -  6a2b2<rU2 


+  b4  0"  2 ]  •  [a4  -  6a2b2  +b4]  +  [-4uv +  4u3v^p 

v^c r  2  -  4uv 


2.2  .4 


-  6a 

2^2  ,  3  -  ^  2 


-  4uv3^2  -  4uv602^CT2  +  4u3v€CT2  -  4uv3^^  2 


+  C0 


2tfp«  +  a^afe*-  2o^ffp: 


+  2  602v2  tfn2  + 


u4  Up  2  _  6u2v2  #p  2  +  V4  Sp  2  +  Cl>2  ^0"  2  +  CO  2  *^2  5cr  2 
-  >?oo2u2  if2(T2  +  2oj2v2  ffO-2  +  U4  er 0~ 2  -  6u2v2  $Q-2 


+  v43"0"2]’[u4  “  6u^v^  +  v4]  +  [4u'3v  -  4uvJ] 


2  2  4 


t^e2t°2  + 

6^V,2  +  +  ^o-2^  +  ^cr26o2V 


2  2  2 


2  ,  .2  2 


2  .2  2 
.v 


2 ,2  2 


2«o-2w2u2  +  z^u>y  *  trV  -  e^o-'u 


.2.  .2  2 


2  4 


2  2  2 
v 


2  3  „  x  n2  3 

v  +  4  0  D  uv  + 


+  ^0“  2v^  +  4  p2uv6<i2  -  4  US  p  2u' 

4  UCT^uvcOg  “  4^CT2u3v  +  4?JG"2uv3]  +  [4a3b  -  4ab3]  * 

2,  .4  ^  ^  2  ,  „2  2 


t  +  e^aWg^  -  2p^a2  b»  * 


2  .  2 
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-  6^2aV  +  ^  V  +  2^2  +  p>tf2^02^2  - 

^f^2^2a2  +  Y  2C^b2  +  p  V  -  6  p  Y2a2b2  + 

p  U  2b^  +  4{J^2602  ab  -  4CJ“^  2a3b  +  4  Q*^  2ab3 
+  4CT3  2602ab  -  4  (rtf  2a3b  +  4C*1$  2ab3] 

Computations  involving  Equation  (53)  were  also  carried  out 
on  a  high  speed  computer. 
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3 .  RESULTS 


The  results  of  this  paper  are  presented  in  both  tabular 
and  graphical  forms.  For  all  cases  considered,  and  K  are 
set  equal  to  unity. 

Table  1  lists  values  of  the  mass  ratio  M  of  0.02,  0.05, 

0.1,  and  0.2;  typical  values  of  the  loss  factor  of  0.2. 

0.5,  and  1.0;  and  values  of  the  natural  frequency  to  ^  of  the 

primary  system  of  10,  20,  30,  40,  and  50  radians  per  second. 

For  each  combination  of  the  above  parameters,  there  are  listed 

two  values  of  the  natural  frequency  60  2  °f  the  damper  unit. 

The  first  column  of  60  o  represents  the  optimum  CO  Q  for 

Z  2 

minimum  mean  square  response  when  the  system  is  excited  by  a 
white  noise  at  n^.  The  second  column  of  60  2  represents  the 
optimum  LO  2  f°r  minimum  mean  square  response  when  the  foun¬ 
dation  is  moved  in  a  random  fashion. 

Table  2  lists  values  of  the  mass  ratio  M  of  0.05,  0.1, 
and  0.2;  values  of  the  loss  factor  of  0.2  and  0.5;  and 
values  of  the  natural  frequency  60  of  the  primary  system 
of  10,  20,  30,  40,  and  50  radians  per  second.  The  first 
column  of  60  2  represents  the  optimum  60  2  for  the  equal 
peak  criterion  when  the  excitation  is  applied  at  The 

second  column  of  60  2  represents  the  optimum  £0  2  for  the 
equal  peak  criterion  when  the  excitation  is  the  motion  of 
the  foundation.  Data  for  the  equal  peak  criterion  was  obtained 
for  only  a  limited  number  of  combinations  of  the  parameters. 

For  the  combinations  of  parameters  listed  in  Table  1  that  are 
not  listed  in  Table  2,  the  absolute  value  of  the  frequency 
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response  function  exhibited  essentially  only  one  peak.  The 
cases  where  all  four  sets  of  optimum  U)2  were  obtained  are 
plotted  in  Figures  4,  5,  6,  and  7. 

To  further  contrast  the  different  results  obtained  by 
use  of  two  different  tuning  criteria,  in  Figures  b,  9  and  10 
are  plotted  the  square  of  the  absolute  value  of  the  fre¬ 
quency  response  function  versus  60  for  a  mass  ratio  M  =  0.05, 
a  loss  factor‘d  =  0.2,  and  a  natural  frequency  of  the  primary 
system  60  ^  =  40  radians  per  second.  Figure  8  represents  the 
case  of  60  2  =  38.4  radians  per  second,  which  is  the  value  of 
60  g  that  satisfies  the  equal  peak  criterion.  Figure  9 
represents  the  case  of  60  2  =  38.9  radians  per  second,  which 
is  the  value  of  60  2  that  satisfies  the  minimum  mean  square 
response  criterion  when  the  system  is  excited  by  a  white 
noise  at  m1.  This  value  of  GO  2  causes  the  first  peak  to  be 
of  greater  amplitude  than  the  second  peak.  Figure  10  repre¬ 
sents  the  case  of  60  2  =  38.0  radians  per  second.  This 
value  of  60  2  is  less  than  both  the  frequency  values  corre¬ 
sponding  to  the  two  tuning  criteria  stated  above,  and  it 
causes  the  second  peak  to  be  of  greater  amplitude  than  the 
first  peak. 

Table  3  lists  values  of  loss  factor  of  0.07,  0.10, 
0.20,  0.50,  and  1.0  in  combination  with  a  single  mass  ratio 
M  =  0.02  and  a  single  natural  frequency  of  the  primary  system 
601  =  50  radians  per  second.  For  each  value,  nine  values 
of  60  are  listed  together  with  their  respective  values  of 
mean  square  response  for  the  case  of  white  noise  excitation 
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2 

at  m^.  The  units  on  E[X  (t)]  would  depend  on  the  units  used 

for  the  masses  in  the  system.  If  m^  and  m ^  are  specified  by 

2  2 

the  units  pound-second  per  inch,  then  E[X  (t)]  would  be 

o 

expressed  in  terms  of  inch  .  The  data  from  Table  3  is  plotted 
in  Figure  11. 

Table  4  contains  the  same  combinations  of  parameters  as 
Table  3,  except  that  Table  4  refers  to  the  case  of  random 
motion  of  the  foundation.  The  data  of  Table  4  is  plotted  in 
Figure  12. 

Table  5  lists  values  of  the  mass  ratio  M  of  0,02,  0,05, 
0.10,  and  0,20  in  combination  with  values  of  loss  factor  ^ 
of  0.04,  0.07,  0.10,  0.20,  0.50,  and  1.00,  and  a  single 
CO  ^  s  10  radians  per  second.  For  each  mass  ratio  and  the 
given  00  there  exists  an  optimum  loss  factor  giving  rise 
to  a  minimum  value  of  mean  square  response  for  each  case  of 
white  noise  excitation.  The  data  for  the  case  where  white 
noise  excitation  is  applied  at  m1  is  plotted  in  Figure  13, 
and  the  data  for  the  case  where  white  noise  excitation  is 
applied  at  the  foundation  is  plotted  in  Figure  14. 

Table  6  is  identical  to  that  of  Table  5  except  that  the 
natural  frequency  00x  of  the  primary  system  is  50  radians 
per  second.  The  data  from  Table  6  is  plotted  in  Figures  15 
and  16. 
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4.  CONCLUSIONS 


Figures  5,  6,  and  7  show  that  when  the  excitation  is 
applied  at  the  two  criteria  for  optimization  give  rise  to 

optimum  values  of  CO  2  that  are  not  identical,  but  very  close 
to  each  other.  When  the  excitation  is  the  motion  of  the  foun¬ 
dation,  the  two  criteria  also  give  rise  to  optimum  values  of 
CO g  that  are  very  close  to  each  other. 

It  is  also  seen  that  the  plots  of  to  1  versus  optimum 
are  essentially  linear.  By  comparing  Figures  6  and  7, 
it  is  found  that  for  a  given  mass  ratio,  the  slope  of  the 
plots  decreases  as  loss  factor  decreases.  This  implies  that 
the  greater  the  loss  factor  the  softer  the  damper  unit  must 
be  to  satisfy  either  optimum  criterion. 

Figures  8,  9,  and  10  show  that  the  second  peak  is  higher 
than  the  first  peak  for  a  value  of  60  2  that  is  smaller  than 
the  one  giving  rise  to  equal  peaks.  The  reverse  is  true  If 
60 ^  is  larger  than  that  which  gives  rise  to  equal  peaks.  If 
60  g  is  considerably  smaller  than  or  considerably  larger  than 
the  optimum  CO  essentially  only  one  peak  will  appear  in 
the  absolute  value  of  the  frequency  response  function. 

As  mentioned  previously  in  this  section,  Figures  4,  5, 

6,  and  7  show  that  the  two  criteria  for  optimization  are  not 
exactly  identical.  Figures  11  and  12,  however,  show  that  the 
plots  representing  mean  square  response  versus  60  ^  are  Quite 
flat  in  the  region  of  their  respective  minima.  It  can  be  seen 
that  for  a  given  value  of  mass  ratio  and  a  given  value  of  60  , 

the  greater  the  loss  factor,  the  more  slowly  the  mean  square 
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response  varies  in  the  region  of  its  minimum.  For  the  case 
of  mass  ratio  M  =  0.2,  loss  factor  ^  =  0.2,  and  60  ^  =  50 
radians  per  second,  the  equal  peak  condition  is  accompanied 
by  a  mean  square  response  which  is  only  two  per  cent  larger 
than  the  minimum  mean  square  response. 

Figures  11  and  12  show  that  as  the  value  of  60  2  deviates 
from  the  optimum  value  of  iOgj  the  mean  square  response 
increases.  For  the  limiting  cases,  that  is,  60  2  =  0  and 
60  j  i  the  mean  square  response  is  unbounded.  The  case 

of  to  g  c  ®  corresponds  to  the  case  where  the  damper  unit  is 
not  attached  to  the  primary  system.  Since  the  primary  system 
is  of  a  single  degree-of-freedom  in  which  no  damping  is  pres¬ 
ent,  the  frequency  response  function  becomes  unbounded  at 
CO  e  60  For  this  case  the  area  under  the  squared  frequency 
response  function  curve  is  also  infinite.  For  the  case  of 
CO  2  =  oe  ,  the  damper  spring  is  infinitely  stiff  so  that  the 
viscoelastic  link  becomes  a  rigid  link,  and  again  the  system 
degenerates  into  a  single  degree-of-freedom  system.  The  mass 
of  this  system  is  made  up  of  plus  m2  Joined  together  by 
the  rigid  link.  Similar  to  the  case  of  60  2  =  0,  the  squared 
frequency  response  function  becomes  unbounded  at  some  frequency 
and  results  in  an  infinte  area  under  the  curve. 

Referring  to  Figures  13,  14,  15,  and  16  it  can  be  seen 
that  for  each  combination  of  mass  ratio  and  60  ^  there  exists 
a  value  of  the  loss  factor  that  will  give  rise  to  a  minimum 
mean  square  response.  For  a  given  value  of  CO  the  optimum 
loss  factor  decreases  as  the  mass  ratio  decreases.  The 
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optimum  loss  factor  for  a  given  mass  ratio  does  not  appear 
to  be  a  function  of  60  1  as  can  be  seen  by  comparing  Figure 
12  with  Figure  14  and  Figure  13  with  Figure  15.  Data  for 
the  case  of  60  ^  =  30  radians  per  second  also  verifies  this. 
However,  the  value  of  the  minimum  mean  square  response 
corresponding  to  the  optimum  loss  factor  does  appear  to  be 
a  function  of  CO  This  can  be  seen  by  again  comparing 
Figures  12  and  14  and  Figures  13  and  15. 

In  short,  it  can  be  said  that  for  a  combination  of  loss 
factor,  mass  ratio,  and  60  such  that  the  equal  peak  condi 
tion  can  be  satisfied,  the  optimum  C02  corresponding  to 
equal  peaks  is  almost  equal  to  that  which  satisfies  the  mini 
mum  mean  square  response  criterion.  However,  it  should  be 
noted  that  there  are  combinations  of  the  above  parameters 
that  allow  the  system  to  be  optimized  according  to  minimum 
mean  square  response  but  do  not  allow  the  system  to  be 
optimized  according  to  the  equal  peak  criterion. 
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TABLE  1 


Optimum  Values  of  OJ2  for  Minimum  Mean  Square  Response 


M  V  ^1 


0.02  0.2  10 
20 
30 
40 
50 

0.5  10 

20 
30 
40 
50 

1.0  10 
20 
30 
40 
50 

0.05  0.2  10 

20 
30 
40 
50 

0.5  10 

20 
30 
40 
50 

1.0  10 
20 
30 
40 
50 


2  f°r  white  noise 
at  m^ 

9.80 

19.6 

29.4 
39.3 

49.2 

9.36 

18.6 

28.3 

37.7 

47.1 

8.37 

16.5 
25.0 

33.8 

43.5 
9.74 

19.5 

29.2 

38.9 

47.8 

9.22 

18.5 

27.9 

37. 1 

46.8 

8.22 

16.3 

24.8 

33.1 

41.4 


^2  f°r  white  noise 
at  foundation 

9.95 

19.9 

29.8 

39.7 

49.5 
9.49 

19.0 

28.6 

38.2 

47.8 

8.39 

16.6 

25.2 

33.9 

44.2 

10.0 

19.9 

29.9 
39.75 

49.6 
9.60 

19. 1 

28.7 

38.2 

47.9 

8.40 

16.7 

25.3 

33.6 
42.5 
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TABLE  1  (continued) 


Optimum  values  of  CO  „  for  Minimum  Mean  Square  Response 


M  ij  60  j 


0.1  0.2  10 
20 
30 
40 
50 

0. 5  10 

20 
30 
40 
50 

1.0  10 
20 
30 
40 
50 

0.2  0.2  10 
20 
30 
40 
50 

0.5  10 

20 
30 
40 
50 

1.0  10 
20 
30 
40 
50 


60g  for  white  noise 
at 

9.65 

19.2 

28.7 

38.4 

47. 1 
9.05 

18.2 

27.4 

36.2 

45.6 
8.12 

16.3 

24.2 
32,  4 
41.0 

9.28 

18.6 
27. 9 

37 . 1 

46.5 
8.90 

17.7 

26.6 

35.4 

43.7 

8.00 

15.8 

23.9 

31.1 

38.8 


CO  for  white  noise 
at  foundation 

10.05 

20.2 

30.1 

40.1 

50.1 
9.60 

19.2 

28.8 

38  4 
48.0 

8.52 

16.9 

25.7 

33.9 

42.7 

10.2 

20.5 

30.8 

41,2 

51 . 2 
9.70 

19.6 
29  5 

39  3 
49.0 

8.66 

17.3 

25.8 

34.2 

43.4 
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TABLE  2 


Optimum  Values  of  CO^  for  Equal  Peak  Criterion 


1 

60 

1 

60u  for  sinusoidal 
input  at 

<«0  for  sinusoidal 
input  at 
foundation 

0.2 

10 

9.60 

9.90 

20 

19.2 

19.8 

30 

28.7 

29.6 

40 

38.4 

39.5 

50 

48.0 

49.4 

0.2 

10 

9.40 

10,0 

20 

18.8 

19  9 

30 

28.1 

29.7 

40 

37.6 

39.6 

50 

46.9 

49.5 

0.2 

10 

9.  10 

9.90 

20 

17.9 

19.9 

30 

26.7 

29.8 

40 

35.9 

39.7 

50 

44.9 

49.6 

0.5 

10 

8.50 

9.  50 

20 

16.9 

18.7 

30 

25.4 

27.9 

40 

33.9 

37.3 

50 

42.6 

46.9 
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TABLE  3 


Mean  Square  Response  for  the  Case  of  White  Noise 

Excitation  at  m. 


0.02 


1 


0.07 


50 


0.10 


50 


0.20 


50 


0.50 


50 


1.00 


50 


<0 

2 

37.4 

40.0 

42.45 

44.75 

46.9 
49.0 
51.0 

52.9 

54.8 
37.4 
40.0 
42.45 
44.75 

46.9 
49.0 
51.0 

52.9 

54.8 
37.4 
40.0 
42.45 
44.75 

46.9 
49.0 
51.0 

52.9 

54.8 
37.4 
40.0 
42.45 
44.75 

46.9 
49.0 
51.0 

52.9 

54.8 
37.4 
40.0 
42.45 
44.75 

46.9 
49.0 
51.0 

52.9 
54.8 


E[X2(t)]  x  103 


4.97 
3.92 
3.00 
1.01 
0.68 
0.46 
0.51 
0,71 
1.47 
3.88 
2.85 
1.99 
0.84 
0.53 
0.38 
0.42 
0.60 
0.97 
1.58 
1.43 
0.98 
0.61 
0.43 
0.37 
0.41 
0.51 
0.67 
1.06 
0 . 93 
0.76 
0.67 
0.65 
0.67 
0.72 
0.81 
0.93 
1.15 
1.11 


07 

08 

11 

16 

,23 

,30 

.38 
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TABLE  4 


Mean  Square  Response  for  the  Case  of  White  Noise 
Excitation  at  the  Foundation 


M 

0.02 


1 

0.07 


0.10 


0.20 


0.50 


1.00 


60 

1 

<°2 

E[Z2(t)]  x 

50 

37.4 

40.0 

42.45 

2.03 

44.75 

0.71 

46.9 

0.48 

49.0 

0.30 

51.0 

0  298 

52.9 

0  394 

54.8 

0.797 

50 

37.4 

40.0 

42.45 

1.34 

44 . 75 

0.58 

46.9 

0.365 

49.0 

0.246 

51.0 

0.245 

52.9 

0.334 

54.8 

0.533 

50 

37.4 

1.03 

40.0 

0.95 

42.45 

0.65 

44.75 

0.40 

46.9 

0.28 

49.0 

0.23 

51.0 

0.24 

52.9 

0 . 29 

54.8 

0.38 

50 

37.4 

0.76 

40.0 

0 , 59 

42.45 

0.48 

44.75 

0.42 

46.9 

0.40 

49.0 

0.41 

51.0 

0.43 

52.9 

0.48 

54.8 

0. 55 

50 

37.4 

0.76 

40.0 

0.68 

42.45 

0.654 

44.75 

0.652 

46.9 

0.67 

49.0 

0.70 

51.0 

0.73 

52.9 

0.78 

54.8 

0.82 
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TABLE  5 


Minimum 

M 

0.02 

0.05 

0.10 

0.20 


Mean  Square  Response  for  Both  Cases  of  White  Noise 
Excitation 


T 

“i 

E [ X2 ( t ) ]  x  101 
excitation  at 

E[Z2  (t )  ]  x 
excitation 
foundation 

0.04 

10 

0.916 

0.921 

0.07 

0.56 

0.58 

0.10 

0.47 

0.488 

0.20 

0.473 

0.485 

0.50 

0.816 

0  818 

1.00 

1.35 

1.34 

0.04 

10 

0.816 

0,822 

0.07 

0.501 

0.514 

0.10 

0.385 

0.396 

0.20 

0,283 

0.298 

0.50 

0.369 

0.372 

1.00 

0.568 

0  562 

0.04 

10 

0.785 

0  856 

0.07 

0.467 

0.499 

0.10 

0.343 

0  369 

0.20 

0.220 

0  236 

0.50 

0.219 

0  223 

1.00 

0.309 

0.301 

0.20 

10 

0.186 

0  207 

0.50 

0.  144 

0.  148 

1.00 

0.178 

0.  169 
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TABLE  6 


Minimum  Mean  Square  Response  for  Both  Cases  of  White  Noise 

Excitation 


M 

"*L 

^1 

E[X2(t)]  x  103 

E[Z2 (t) ]  x 

excitation  at  m. 

excitation 

1 

foundation 

0.02 

0.04 

50 

0.63 

0.39 

0.07 

0.44 

0.27 

0.10 

0.38 

0.23 

0.20 

0.37 

0.23 

0.50 

0.65 

0,40 

1.00 

1.07 

0.65 

0.05 

0.04 

50 

0.61 

0.41 

0.07 

0.38 

0  24 

0.10 

0.29 

0.186 

0.20 

0.22 

0. 137 

0.50 

0.29 

0  175 

1.00 

0.45 

0,269 

0. 10 

0.04 

50 

0.60 

0.39 

0.07 

0.37 

0.23 

0.10 

0.27 

0  172 

0.20 

0.  172 

0,  105 

0.50 

0.174 

0.  09 

1.00 

0.  244 

0, 138 

0.20 

0.20 

50 

0.  146 

0. 087 

0.  50 

0. 112 

0.06 

1.00 

0.14 

0.072 
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Figure  1.  Random  excitation  at 
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Figure  2.  Random  motion  of  foundation 
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Figure  3.  Location  of  poles  and  path  used  in 
contour  integration 
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60^  =  natural  frequency  of  primary  system 
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Figure  6,  Plot  of  40^  versus  optimum  £0* 


34 


50 

40 

30 

20 

10 

0 


m  ^<3  o  O  •  3  3 


=  0.2 

=  0.5 

-  random  excitation  at  m, 

-  random  excitation  at  foundation 
"  sinusoidal  excitation  at  m.. 

-  sinusoidal  excitation  at  foundation 

1  =  natural  frequency  of  primary  system 

2  =  natural  frequency  of  damper  unit 


— i - 1 - 1 - «— 

10  20  30  40 

6^2 (rad/sec) 


Figure  7. 


Plot  of 


versus  optimum  60^ 


— i- 

50 


35 


.05 


c 

o 


o 

CJ  X 

O 

in  ^ 

X 

2  3 

3  m 

D 

tn 

u 

Q) 

> 

CNJ 


3 

s/ 

PC 


O 


O 

r-t 

0< 

oo 

<u 

a 

be 


36 


Sh 


o 

to 


o 

in 


o 


o 

cn 


o 

CM 


o 


gOI  X  z  )(co)H | 


37 


^  (rad/sec) 

Figure  9.  Plot  of  |h(^)|  2  versus  ti)  (excitation  at 
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^  (rad/sec) 

versus  60  (excitation  at 


E  [  X  ■(  t )  ] 


35 
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Figure  11.  Plot  of  mean  square  response  E^X  (t)] 
versus  natural  frequency  of  the  damper 
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E[Zz(t)]x  10 
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Figure  12.  Plot  of  mean  square  response  E[Z  tt)) 
versus  natural  frequency  of  damper 
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Figure  14,  Plot  of  minimum  mean  square  response  E[Z^(t)] 
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Figure  15,  Plot  of  minimum  mean  square  response  E[X  (t)] 
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